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ABSTRACT 

We present a catalogue of 179 hyperluminous infrared galaxies (HLIRGs) 
from the Imperial IRAS-FSS Redshift (IIFSCz) Catalogue. Of the 92 with 
detections in at least two far infrared bands, 62 are dominated by an M82- 
like starburst, 22 by an Arp220-like starburst and 8 by an AGN dust torus. 
On the basis of previous gravitational lensing studies and an examination 
of HST archive images for a further 5 objects, we estimate the fraction of 
HLIRGs that are significantly lensed to be 10-30%. 

We show simple infrared template fits to the SEDs of 23 HLIRGs 
with spectroscopic redshifts and at least 5 photometric bands. Most can 
be fitted with a combination of two simple templates: an AGN dust torus 
and an M82-like starburst. In the optical, 17 of the objects are fitted with 
QSO templates, several with quite strong extinction. There are 5 objects 
fitted with galaxy templates in the optical, two of which show evidence for 
AGN dust tori and so presumably contain Type 2 (edge-on) QSOs. The 
remaining object is fitted with a galaxy template in the optical, but is of 
such high luminosity that this classification would be plausible only if there 
were very strong lensing. 20 of the 23 objects (87%) show evidence of an 
AGN either from the optical continuum or from the signature of an AGN 
dust torus, but the starburst component is the dominant contribution to 
bolometric luminosity in 14 out of 23 objects (61 %). The implied star- 
formation rates, even after correcting for lensing magnification, are in 
excess of 1000 M^yr^ 1 . 

We use infrared template-fitting models to predict fluxes for all 
HLIRGs at submillimetre wavelengths, and show predictions at 350 and 
850 fim. Most would have 850 fim fluxes brighter than 5 mJy so should 
be easily detectable with current submillimetre telescopes. At least 15 % 
should be detectable in the Planck all-sky survey at 350 /im and all Planck 
all-sky survey sources with z < 0.9 should be IIFSCz sources. 

From the luminosity-volume test we find that HLIRGs show strong 
evolution. A simple exponential luminosity evolution applied to all 
HLIRGs would be consistent with the luminosity functions found in red- 
shift bins 0.3-0.5, 0.5-1 and 1-2. The evolution-corrected luminosity func- 
tion flattens towards higher luminosities perhaps indicating a different 
physical mechanism is at work compared to lower luminosity starbursts. In 
principle this could be gravitational lensing though previous searches with 
HST have, perhaps surprisingly, not shown lensing to be widely prevalent 
in HLIRGs. 

Key words: infrared: galaxies - galaxies: evolution - stardormation - 
galaxies: starburst - cosmology: observations 
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1 INTRODUCTION 

Rowan-Robinson (2000, hereafter RR2000) defined 
hyperluminous infrared galaxies to be those with 
bolometric infrared luminosities L; r > lO 13 !/©. He 
discussed a sample of 39 such galaxies found either by 
follow-up of the IRAS survey or through submillime- 
tre observations, and modelled the infrared and sub- 
millimetre spectral energy distributions (SEDs) for 
those that were well-observed. Submillimetre obser- 
vations and more detailed SED models for 13 hyper- 
luminous infrared galaxies (HLIRGs) were reported 
by Farrah et al (2002a) and ISO data for 4 HLIRGs 
was discussed by Verma et al (2002). A more de- 
tailed SED model for the prototypical HLIRG IRAS 
F10214+4724 was give by Efstathiou (2006). Recently 
Ruiz et al (2010) have modelled the SEDs of 13 hy- 
perluminous galaxies which are X-ray sources. Many 
of these galaxies show evidence of AGN dust torus 
emission at rest-frame wavelengths 3-30 pm, but the 
submillimetre emission is almost always due to star- 
formation, an interpretation that is supported by the 
large gas-masses found in these galaxies (see RR2000 
for a summary). 7 of the 39 objects in RR2000 were 
already known to be gravitationally lensed. A further 
9 hyperluminous infrared galaxies were imaged with 
HST by Farrah et al (2002b), but none were found 
to be lensed. Even after correction for magnification 
by gravitational lensing, most of these galaxies imply 
star formation rates > lOOOMQi/r -1 . 

Since many of the hyperluminous infrared galax- 
ies in RR2000 were discovered through their submil- 
limetre radiation, it is clearly of interest to consider 
whether all or most submillimetre galaxies are hy- 
perluminous. Many discussions of the nature of 'sub- 
millimetre galaxies' assume the latter to be the case. 
However there is growing evidence that the submil- 
limetre population is quite diverse. In their models 
for the SEDs of a complete sample of 850 /im sources 
from the SHADES survey, Clements et al (2008) find 
a diverse range of SED types, including M82 star- 
bursts, 'cirrus' (ie quiescent) galaxies, and AGN dust 
tori, as well as Arp220-type starbursts. Efstathiou 
and Rowan-Robinson (2003) have pointed out that 
the SEDs of as many as half of submillimetre galax- 
ies could be interpreted in terms of a cirrus model. 
In modelling source-counts from 8-1100 /im, and the 
far infrared and submillimetre background radiation, 
Rowan-Robinson (2008) concludes that at S(850) = 8 
mjy there are approximately equal numbers of M82 
and Arp220 starbursts. However cirrus galaxies, in 
these models, dominate the counts only at S(850) > 
100 mjy and at S(850) < 1 mjy. 

The creation of the Imperial IRAS FSC Red- 
shift Catalogue (IIFSCz, Wang and Rowan-Robinson 
2009a) gives us the opportunity to define an almost 
all-sky (|fc| > 20°) catalogue of hyperluminous in- 
frared galaxies, sampled at 60 /im. This allows us to 
study the evolution and luminosity function of this 



population and to predict what will be seen in sub- 
millimetre surveys to be undertaken by Planck and 
Herschel. 

A spatially flat cosmological model with A = 0.7, 
/io=0.72 has been used throughout. 



2 THE CATALOGUE 

Our catalogue of hyperluminous IRAS galax- 
ies, with Li r > lO 13 Z/0, derived from IIF- 
SCz, contains 179 galaxies. The catalogue, together 
with the full IIFSCz Catalogue, is available at 
http://astro.imperial.ac.uk/ ~mrr/fss/hyper. 

Of the 25 hyperluminous galaxies in Tables 1 and 
2 of RR2000, found from 60 or 850 /im surveys, 5 were 
too weak at 60 fim to be detected in the FSS, and 2 lie 
at |b| < 20 . Of the remaining 18, 14 are in the present 
catalogue, 3 lie just below the 10 13 L Q limit, and one 
has had its redshift revised to a lower value. On the 
basis of the 12 hyperluminous galaxies detected from 
unbiased IRAS surveys, Rowan-Robinson (2000) esti- 
mated that there would be 100-200 such galaxies with 
S(60) > 200 mjy over the whole sky . 

Figure 1 shows a plot of Li r versus redshift for 
IIFSCz galaxies, colour-coded by the dominant in- 
frared template type. Note that in fitting infrared 
templates we specifically exclude the possibility of an 
HLIRG having a cirrus template. Of the 92 galaxies 
with detections in at least two infrared bands, the 
infrared luminosities of 62 are dominated by an M82- 
like starburst, 22 by an Arp 220-like starburst, and 8 
by an AGN dust torus. Where only one infrared band 
is available we have estimated the luminosity using 
an M82 template. 

We have spectroscopic redshifts from the liter- 
ature for 58 (32%) HLIRGs, photometric redshifts 
for the remainder. We have spectroscopic redshifts 
for 18 out of 46 (39%) of galaxies with infrared lu- 
minosities, Li r > lO 14 Z/0. Four hyperluminous in- 
frared galaxies from this catalogue are known to be 
gravitationally lensed: F08105+2554=HS0810+2554 
(Reimers et al 2002), F08279+5255 (Downes et al 
1998), F10214+4724 (Eisenhardt et al 1996) and 
F14132+1144=H1413+117 (Yun et al 1997). The first 
of these does not have a lens model, but the other 
three remain in the hyperluminous category when cor- 
rected for their estimated far infrared lensing magni- 
fications of 14, 10 and 10, respectively (see Table 1). 
We have adopted a far infrared lensing magnification 
of 10 for F08105+2554 by analogy with the very sim- 
ilar system PG1115+080 (Reimers 2002). A further 
six hyperluminous galaxies were found by Farrah et al 
(2002) not to be lensed (F00235+1024, F10026+4949, 
F12509+3122, F14218+3845, F15307+3252, 
F16348+7037=PG1634+706). We examined the HST 
archive images for a further five hyperluminous 
galaxies ( F01175-2025, F09105+4108, F14165+0642, 
F09105+4108, F18216+6419) and found no evidence 
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Figure 1. Bolomctric infrared luminosity versus redshift for IIFSCz galaxies. Colour coding indicates dominant contribu- 
tion to bolometric infrared luminosity. Black: cirrus (quiescent); red: M82 starburst; green: Arp220 starburst; blue: AGN 
dust torus; cyan: M82 starburst fitted to 60 /an only. Known gravitationally lensed systems are labelled L. 



of lensing, though F09105+4108 and F09105+4108 
have highly disturbed images indicating a major 
merger. Thus in the sample for which we have infor- 
mation on lensing, the fraction which are significantly 
lensed is 4/15 (27%). On the basis of the relatively 
unbiassed search for lensing in HLIRGs by Farrah et 
al (2002), the fraction of HLIRGs which are signifi- 
cantly lensed may be no greater than 10%. However 
this is clearly a very uncertain estimate and it would 
be worthwhile searching for new lenses in this sam- 
ple, especially amongst the galaxies with infrared lu- 
minosities, Li r > lO 14 Z/0. The lensing rate might be 
expected to be higher for high redshift (z>2) galaxies. 
There are 20 such galaxies in our catalogue and 3 are 
already known to be lensed. 



3 SPECTRAL ENERGY DISTRIBUTIONS 
AND PREDICTED SUBMILLIMETRE 
FLUXES 

Detailed modelling of SEDs of hyperluminous galaxies 
was presented by Rowan-Robinson (2000) and Farrah 
et al (2002b), and the latter presented some new sub- 
millimetre observations for 13 hyperluminous galax- 
ies. Efstathiou (2006) has modelled the prototypi- 
cal hyperluminous galaxy F10214+4724, for which 
Spitzer IRS data was obtained by Teplitz et al (2006) . 
Figures 2-4 show SEDs for 23 hyperluminous galax- 
ies, which have spectroscopic redshifts and photomet- 
ric data in at least 5 bands, including at least two 
infrared bands. Data have been compiled from the 
literature. Table 1 gives physical parameters for the 
23 sources. The columns are: name, redshift, lensing 
magnification, luminosities in starburst, dust torus, 
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Figure 2. SEDs for 8 hypcrluminous galaxies with spectroscopic rcdshifts. The template fits used in the infrared are M82 
star burst (dotted curve), AGN dust torus (long-dashed curve). Parameters for fits given in Table 1. 



and optical galaxy/QSO components, optical tem- 
plate type, visual extinction, star-formation rate. Lu- 
minosities in this table have been corrected for the 
estimated lensing magnification. 12 have strong AGN 
dust tori, while 9 show no evidence for a dust torus, 
though this may be partly because of a lack of obser- 
vational data at rest-frame wavelengths 3-30 /jm . In 
fact the infrared SEDs of all 23 objects can be fitted 
with mixtures of just two components, an AGN dust 
torus and an M82 starburst. The model fits in Figs 2-4 
broadly agree with the automatic template fits just to 
IRAS data: the relative luminosities of starburst and 
dust torus are slightly modified in some cases when 
fitting to the additional far infrared and submillime- 
tre photometry here. No galaxies with Arp 220 fits 
in the IIFSCz Catalogue had a sufficient number of 
photometric bands to be plotted in these figures. 

In the optical, 17 of the objects are fitted with 



QSO templates, three with strong extinction (Ay > 
0.5). There are 5 objects fitted with galaxy templates 
in the optical, two of which show evidence for AGN 
dust tori and so presumably contain Type 2 (edge- 
on) QSOs. The remaining object (F08105+2554) is 
best fitted with a galaxy template in the optical, but 
is of such high optical luminosity that this classifica- 
tion would be plausible only if there were very strong 
lensing. As this is a known quadruple lensed system 
(Reimers et al 2002), this is a possibility (our estimate 
of the lensing magnification is 10 - see section 2): a 
QSO fit is also shown, which fits the optical data but 
not the JHK data. For F14218+3845 there is a dis- 
crepancy between the IRAS fluxes and the 90 fim flux 
measured with ISO. Our preferred model is the more 
conservative fit to the ISO data. 20 of the 23 objects 
(87%) show evidence of an AGN either from the opti- 
cal continuum or from the signature of an AGN dust 
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Figure 3. SEDs for 8 hyperluminous galaxies. Parameters for fits given in Table 1. The IRS data for F10214+4724 are 
indicated in red and the Efstathiou (2006) torus model for this object is shown as a short-dashed locus. 



torus, but the starburst component is the dominant 
contribution to bolometric luminosity in 14 out of 23 
objects (61 %). Even if the QSO is not the dominant 
source of the extreme infrared luminosities, there is 
clearly a strong case for a link between the very high 
rates of star-formation implied by these galaxies and 
the growth of massive black holes (cf Sanders et al 
1988, Kauffmann and Haehnelt 2000, Vielleux et al 
2009). However it should be borne in mind that re- 
stricting the sample of Fig 2 and Table 1 to objects 
with spectroscopic redshifts inevitably biases the sam- 
ple towards QSOs. 

Our simple dust torus template does not fit the 
data for F10214+4724, which Efstathiou (2006) sug- 
gests is a highly edge-on system. We have used the 
Efstathiou (2006) edge-on dust-torus model in our 
fit to F10214+4724 (Fig 3). The J,H,K data for 
F15415+1633 are brighter than our model predicts, 



and this is true also for F08105+2554 if a QSO model 
is preferred for the optical SED. Three strong radio 
sources which have very flat SEDs and are strongly 
variable, OJ287, 3C345 and 3C446, have not been in- 
cluded in the SED plots. 

For QSOs the ratio of torus luminosity to optical 
luminosity, L t0 r/L op t, can be interpreted as the dust 
covering factor (Rowan- Robinson et al 2009). In this 
sample there are several objects with QSO-like optical 
SEDs for which L t0 r is significantly greater than L opt 
(F10026+4949, F10214+4724, F15307+3253) and for 
these we must assume that we are seeing only a small 
fraction of the QSO optical light, probably via scat- 
tering off the dust torus (cf Polletta et al 2006) . 

To test the validity of our template fits to the 
IRAS data we show in Fig 5 composite spectra of 
HLIRGs, detected in at least two far infrared bands, 
which select M82 or Arp 220 templates. The distinc- 
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Figure 4. SEDs for 7 hypcrluminous galaxies. Parameters for fits given in Table 1. 



tion between these two SED types appears to be gen- 
uine. 

Wang and Rowan-Robinson (2009a) fitted four 
infrared templates to the far infrared SEDs of the 
whole IIFSCz Catalogue and used these to predict 
submillimetre fluxes at a range of wavelengths and 
we have used these template fits to predict 350 and 
850 pm fluxes, and to investigate evolution. Efs- 
tathiou and Rowan-Robinson (2009) have analyzed 
Spitzer-IRS mid-infrared spectroscopic data for in- 
frared galaxies, in particular the diagnostic diagram of 
Spoon ct al (2006), in terms of sequences of starburst, 
cirrus (quiescent), and AGN dust tori models. They 
conclude that infrared template fitting could be im- 
proved by including additional templates correspond- 
ing to very young and very old starbursts. However 
the shortage of far infrared photometric data (rarely 
more than 2 detected bands) for the present sample 
means this would not be very meaningful at this stage. 



We have used our 4-template fits to predict sub- 
millimetre fluxes for the IIFSCz catalogue. Figure 6 
shows predicted 850 /im flux versus z, with hyper- 
luminous galaxies indicated by open squares. Most of 
the hyperluminous galaxies in our catalogue have pre- 
dicted 850 fim fluxes brighter than 5 mjy and should 
be easily detectable with current submillimetre tele- 
scopes. 

Figure 7 shows predicted 350 mu flux versus z. 
We have indicated the estimated PLANCK all-sky 
survey detection limit. Out to z ~0.9 all sources de- 
tected by Planck would be in the IIFSCz (for |6| > 
20°). At least 15 % of hyperluminous infrared galaxies 
in our catalogue should be detectable by PLANCK. 
Clearly many other luminous infrared galaxies at z > 
1.0 will also be detected by PLANCK. 

Individual flux predictions are very uncertain, 
since they are based on 25-100 /im fluxes and in many 
cases on 60 /im detections only. For example if an M82 
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Figure 5. Composite SEDs for hyperluminous galaxies, detected in at least two far infrared bands, which select either 
M82 or Arp 220 starburst templates. 



starburst is incorrectly classified as an Arp220 type, 
the submillimctrc fluxes and bolometric luminosities 
can be overestimated by an order of magnitude. 



4 EVOLUTION AND LUMINOSITY 
FUNCTIONS 

To investigate whether there is evidence for evolution 
in luminous infrared galaxies we first plot infrared lu- 
minosity versus comoving volume, Fig 8L. If no evo- 
lution was present then in each band of luminosity 
there would be a uniform distribution of objects with 
volume (Rowan- Robinson 1968). Clearly for ULIRGs 
and HLIRGs there is strong evolution. Fig 8R shows 
the corresponding plot corrected for luminosity evo- 
lution with Q = 3.0. The solid curves denote the FSS 
90% completeness limit of S60 = 0.36 Jy and only 



sources brighter than this limit have been included in 
these plots. 

Wang and Rowan-Robinson (2009b) presented 
luminosity functions for IIFSCz galaxies. Here we fo- 
cus on the high luminosity end of the luminosity func- 
tion using the 1/Vmax method. For galaxies without 
an estimate of the best-fitting infrared template, we 
assume a SED of the form S v oc i>~ 2 to derive the K- 
correction term. Figure 9 shows the 60 /im luminosity 
function for IIFSSz galaxies with L 6 o > 10 12 Lq and 
5*60 > 0.36 Jy, divided into three redshift bins, 0.3- 
0.5, 0.5-1.0 and 1.0-2.0. The errors are derived from 
a = Y,i^ ?t ght ' l , where each galaxy is weighted by 

the redshift completeness as a function of 60 micron 
flux density. Clearly the local 60 micron luminosity 
function (Wang & Rowan-Robinson 2009b) indicated 
by the broken curve is not a good description of these 
ultraluminous and hyperluminous objects due to sig- 
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Table 1. Parameters for SED fits 
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nificant evolution. Figure 10 shows the correspond- 
ing luminosity function after correction for luminosity 
evolution of the form 

L(z) = L(0) exp (Q(t - t)/t ) (1) 

with Q = 3.0. This is a higher rate of evolution 
than that deduced for the whole IIFSCz sample, Q = 
1.7 (Wang & Rowan-Robinson 2009b), for z < 0.2. A 
maximum-likelihood approach is also used to estimate 
the evolution of the luminosity function, following 
the procedure described in Wang & Rowan-Robinson 
(2000b). We found a density evolution of the form 
P — 5.1 or a luminosity evolution Q = 2.7 in the 
redshift range 0.3 < z < 2.0. For L 60 < 10 12 - 5 L© the 
luminosity function agrees well with the evolution- 
corrected luminosity function derived at lower lumi- 
nosities, which is of the form derived by Saunders et 
al (1990). The evolution-corrected luminosity func- 
tion then flattens towards higher luminosities, per- 
haps indicating the presence of a new population or 
a different physical mechanism is at work. Gravita- 
tional lensing is one such possible mechanism but we 
have argued above that rather few of the hyperlumi- 
nous galaxies are in fact lensed. The very high star- 
formation rates in HLIRGs may be indicative of major 
merger events between massive galaxies. The analysis 
of star-formation rate versus gas mass deduced from 
CO observations by Rowan- Robinson (2000, Fig 20) 
may indicate a shorter time-scale for the starburst 



(~ 10 7 yrs or less) than in lower luminosity starbursts. 
The space-density of HLIRGs is - 10~ 8 /i? 2 Mpc~ 3 . 



In Wang and Rowan- Robinson (2009a), two 
methods were employed to estimate photometric red- 
shifts. For sources with 2MASS J, H and Ks photom- 
etry, we choose redshifts given by the neural network 
training set method. For sources with SDSS ugriz pho- 
tometry, redshifts derived from the template-fitting 
technique are preferred. In order to investigate the 
impact of photometric redshift errors, we carried out 
a Monte-Carlo analysis by generating random devia- 
tions from the estimated photometric redshifts using 
a Gaussian model with its standard deviation set to 
4.4% in (1 + z). The value 4.4% is chosen because this 
is the worst rms error. We calculated the 1/Vmax lu- 
minosity function of sources with Lqo > W 12 Lq and 
S60 > 0.36 Jy for each of the one hundred Monte 
Carlo simulations created and found that the result- 
ing rms error is several times smaller than the Poisson 
error. The outlier rate is estimated to be only < 0.1%, 
but it might have an impact on our derived luminos- 
ity function, for example, in the highest redshift bin 
(1.0 < 2 < 2.0), where the total number of galax- 
ies in the two luminosity bins is only three and four, 
respectively, of which three and two are photometric 
redshifts. 
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Figure 6. Predicted 850 (im flux versus redshift for IIFSCz galaxies, with hyperluminous galaxies shown as open squares. 
Black: cirrus (quiescent); red: M82 starburst; green: Arp220 starburst; blue: AGN dust torus; cyan: M82 starburst fitted 
to 60 /an only. 



5 CONCLUSIONS 

We present a catalogue of 179 hyperluminous infrared 
galaxies from the Imperial IRAS-FSS Redshift (IIF- 
SCz) Catalogue (Wang and Rowan- Robinson 2009a). 

14 of these were included in the analysis of Rowan- 
Robinson (2000). We use the infrared template-fitting 
models of Wang and Rowan- Robinson (2009a) to pre- 
dict fluxes at 350 and 850 /im. Most would have 850 
/Ltm fluxes brighter than 5 mjy so should be readily de- 
tectable by current submillimetre telescopes. At least 

15 % should be detectable in the Planck all-sky sur- 
vey at 350 fim and all Planck all-sky survey sources 
with z < 0.9 should be IIFSCz sources. 

From the luminosity-volume test we find that 
HLIRGs show strong evolution. A simple exponen- 
tial luminosity evolution applied to all HLIRGs would 
be consistent with the luminosity functions found in 



redshift bins 0.3-0.5, 0.5-1 and 1-2. The evolution- 
corrected luminosity function flattens towards higher 
luminosities perhaps indicating a different physical 
mechanism is at work compared to lower luminos- 
ity starbursts. In principle this could be gravitational 
lensing though previous searches with HST have not 
shown lensing to be widely prevalent in HLIRGs. 

Although for an interesting minority of HLIRGs 
the dominant contribution to the bolometric infrared 
luminosity is due to an AGN dust torus, star forma- 
tion appears to be the dominant factor in the over- 
whelming majority of cases, with star-formation rates 
> 1000 M0j/r _1 being implied. Even if the dura- 
tion of these extreme starbursts is shorter than that 
for less luminous starbursts (Rowan- Robinson 2000), 
they still represent an extremely interesting and dra- 
matic phase of galaxy evolution. 
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Figure 7. Predicted 350 fim flux versus redshift for IIFSCz galaxies, with hyperluminous galaxies shown as open squares. 
Key for coloured symbols as in Fig 2. The horizontal line shows the Planck Blue Book (2005) estimate of the sensitivity of 
the all-sky 350 fim survey (10-cr). 
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Figure 8. Infrared luminosity versus comoving volume for IIFSCz galaxies. LH: no correction for evolution. RH: corrected 
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Figure 9. Upper panel : The 60 fira luminosity function for IIFSSz galaxies with L^q > 1O 12 L0, divided into three redshift 
bins, 0.3-0.5, 0.5-1.0 and 1.0-2.0. The broken curve is the luminosity function calculated by Wang et al (2009) for z < 0.2. 
Lower panel: same, corrected for luminosity evolution of the form eqn (1) with Q = 3.0. 
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